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Abstract 

Quasi-Monte Carlo rules are equal weight quadrature rules defined over the domain 
[0, 1]'*. Here we introduce quasi-Monte Carlo type rules for numerical integration of 
functions defined on M^. These rules are obtained by way of some transformation of 
digital nets such that locally one obtains qMC rules, but at the same time, globally one 
also has the required distribution. We prove that these rules are optimal for numerical 
integration in spaces of bounded fractional variation. The analysis is based on certain 
tilings of the Walsh phase plane. Numerical results demonstrate the efficiency of the 
method. 
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1 Introduction 

Traditionally, quasi-Monte Carlo (qMC) rules are equal weight quadrature formulae defined 
on [0, 1]"^, i.e., 

n=0 
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where P = {xq, . . . , Xjy^i} C [0, 1]* are quadrature points. These rules are used to approxi- 
mate integrals of the form 



see [21 [To] for more information. 

In practice, however, the integrals one needs to approximate are usually over domains 
other than [0, 1]'^, often W. In order for qMC rules to be used in this case, one requires 
a transformation mapping from W to [0,1]** [7]. This step is of great importance, indeed 
different transformations can yield very different results [B]. Furthermore, usually known 
theory does not predict which transformations from W to [0, 1]** will yield the best results [6]. 

To circumvent this problem, we introduce quasi-Monte Carlo type rules for W. Admit- 
tedly, there are still choices to be made by the user of our method. In fact, one could view 
the approach proposed here as a discretized transformation, where the transformed point set 
yields locally qMC rules, but globally has a given distribution. 

To understand what properties a point set locally and globally must satisfy, we analyze 
the worst-case error of numerical integration of functions on W with bounded fractional 
variation of order a and which satisfy a decay condition as one moves away from the origin. 
Roughly speaking, if a = 1 then the functions we consider have square integrable partial 
mixed derivatives up to order one in each coordinate. That is, for a partition V^s of into 
subintervals J = 111=1 have 



The value Vj{f) depends on the function as well as the location J. We assume that for 
J = ni=iK;^i) have Vj(/) — )• as maxi<j<5 min(|aj|, — )• oo with a certain rate 
of decay. Furthermore, also the function / itself must satisfy a certain rate of decay as 
1 1 a? 1 1 — ^ oo (where || ■ || denotes some norm in W). 

For these functions, we now construct quadrature rules for which one obtains an integra- 
tion error bounded above by iV~" (log A^)'^'- (for some Cg^a > depending only on s and 
a), where N denotes the number of quadrature points and s the dimension. The focus in 
this paper is on the asymptotic convergence rate as iV — )■ cxd. The study of tractability [I2l 
is left for future work. 

The design of optimal quadrature rules depends on several properties of the functions 
considered. One is the smoothness of the functions. If the function is concentrated on a 




where 
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bounded region, then the quadrature rule should be able to integrate functions with smooth- 
ness a locally with sufficient accuracy. In other words, one needs enough flexibility to have 
optimal quadrature rules locally. Further, since we can only use a flnite number of quadra- 
ture points iV, for flxed we can only cover a flnite region. Here the assumption of the 
decay of / and the condition on Vj(/) comes in. It tells us where to concentrate our efforts. 
The local property allows us to design quadrature rules which are optimal over subcubes of 
the form ni=i['^«' ^^e global property tells us how important each of those subcubes is 
and how their importance is distributed in the space M*. For example, assume that / and 
Vj(/) decay at least like 11^=1(1 + — )• oo, where x = (xi, . . . ,Xs)- Then the 
most important subcubes are arranged in a hyperbolic cross around the origin. 

Below we show how one can map a digital net deflned on [0, 1]'' such that in each of the 
subcubes, inside this hyperbolic cross region, one obtains a digitally shifted digital net, and 
at the same time, the density of the points in each subcube decreases according to the rate 
of decay of / and Vj(/). Our approach is flexible enough to also be able to handle other 
rates of decay of / and Vj(/). Another way to think about this procedure is the following: 
discretize the transformation from [0, 1]* to such that a digital net in [0, 1]'' is mapped 
to several nets in subcubes of W, such that the subcubes and number of points therein are 
distributed according to a given decay rate. 

The cost of constructing the point set in W is similar to the cost of constructing the 
underlying digital net over a flnite fleld Zf,, assuming that the one- dimensional projections 
have already been deflned. To be more precise, instead of using a mapping [0, 1) as 

is done for digital nets, one uses a look-up table which deflnes a mapping Z]^ i— R. The 
cost of constructing the vectors in Z™ is the same for digital nets in [0, 1)'^ and in W. The 
one-dimensional projections on the other hand need to be carefully designed in advance, as 
we illustrate in Section [3 Since those are designed in advance, we do not count it towards 
the construction cost of an individual point set. 

A numerical test in Matlab, see Section [HI reveals that the computation using the method 
introduced here is considerably faster than using digital nets and the inverse normal cumu- 
lative distribution function (in the numerical example considered here, the proposed method 
is more than ten times faster). The advantage is that in the proposed method one does not 
need to compute the inverse normal cumulative distribution function, which can be time 
consuming. 

Conceivably one could also create the global property by hand: Take a number of digital 
nets Pi, . . . , Pk (for instance from a digital sequence), partition into subcubes Ci, . . . , Ck, 
and put the digital net Pk in Ck ior 1 < k < K. In this case, the global property is simply 
designed by hand (or even adaptively). This method can be become very involved if the 
dimension s is large. Our approach on the other hand builds on a tensor product structure, 
which simplifles matters, especially in high dimensions. To flt, or nearly flt, the integrand 
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to such a structure it may be advantageous to use the method proposed here in conjunction 
with some variance reduction method hke principal component analysis, Brownian bridge 
constructions or similar methods [8]. 

We give an overview of the paper. The analysis of the integration error is based on time- 
frequency analysis. To that end, we introduce tilings of the Walsh phase plane [I5l|16], which 
gives us the flexibility to consider subcubes of different size and adjust them to the rate of 
decay of the function / and the variation Vj(/). This is done in Section [2J In Section [3] we 
analyze the rate of decay of the Walsh coefficients based on the smoothness of /, the rate of 
decay of /, and, roughly speaking, the rate of decay of Vj(/). 

In Section|l]we analyze the integration error: We prove an upper bound on the integration 
error which consists of two parts. One part deals with the remainder, that is, the region 
where there are no quadrature points. For this part one relies on the the rate of decay of 
/ to obtain an upper bound. The other part is concerned with the integration error one 
commits when using a quadrature rule. For this part we use estimates of the integration 
error for qMC rules based on digital nets. 

The construction of the quadrature points is introduced in Section [51 We show how 
one can map a digital net in [0, 1]* into subcubes of M"* such that each subcube contains a 
digitally shifted digital net, where the size of each digitally shifted digital net is according to 
a given distribution. In Section [6] we give upper bounds on the integration error when one 
uses the construction of Section [51 We show that under certain decay rates of / and Vj(/), 
the convergence of the integration error is optimal up to some power of log A^. In Section [7l3l 
we provide three concrete examples of how our bounds can be used in a few situations. 
In Section [HI we represent some numerical results which demonstrate the efficiency of the 
method. 

For the convenience of the reader we introduce some notation used throughout the paper 
in the following subsection. 

Notation 

For s > 1 let S* = {1, . . . , s}. For m C 5 we denote by |m| the number of elements in u. 

The set of complex number is denoted by C, the set of real numbers is denoted by R, 
the set of integers by Z, the set of natural numbers by N, the set of nonnegative integers 
by No, and the finite field of prime order b by Z;,. Further let M.^ = {a; G M : x > 0} and 

= {x G M : X > 0}. For c G C we write c for the complex conjugate of c. For /c, / G Z we 
write k\l if k divides I. 

We always write vectors h = {hi, . . . , hs),j = (ji, . . . ,js), k = (fci, . . . , kg), and so on. 
Further we set = (0, ... , 0) and 1 = (1, . . . , 1) - the dimension of these vectors is apparent 
from the context in which they occur. Further, for u C 5* we set = {hi)i^u and 0) 
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is the vector whose ith coordinate is for i G u and otherwise. For short we write 
foi = (feJi, . . . , b'") and b'^ = {b~^\ . . . , b~^'). For vectors fc, I in and so on, we 

set * / = {kill, . . . , fc^Zs). In particular we write b' -k I = {b^^h, . . . , b^"ls)- Further, for 
/ = {ll, . . . ,ls) we set 

= l^il + ■■■ + |/.|. 

Further we write \x\oo = maxi<j<s and \x\ = {\xi\, . . . , \xs\)- 
For r G Nq we set 

Ur = {i E S : Ti 0}. 
We define the L2 inner product for f,gE L2(M'') as usual by 

{f,9)L2 = / f{x)g{x)dx. 

We say that / is orthogonal to g if (/, g)L2 =0. 

For J C we define the characteristic function by 

I (x) = f ^ if G J, 
otherwise. 

Let P = {xq, . . . , xn-i} C W, a = {\q, . . . , Aat-i} C M. We define the quadrature rule 

N-l 

QpAf) = '^"/(^")- 

n=0 

For A; G N with base b representation k = k,q + Kib + ■ ■ ■ + Ka-ib""^^, we define the vector 
k = (ko, . . . , Km-i)^ G where we set Ka+i = ■ ■ ■ = Km-i = if a < m. 

2 The Walsh model 

For the convenience of the reader we repeat some elementary results concerning the Walsh 
model, see [151 IB] for more detailed information, which is based on Walsh functions, see 

[H [IT]. Let 6 > 2 be an integer, = + 6ki H h t^ab" G No, a; = Xcb" + Xc-ife''"^ H G M 

for some a, c G Nq, be the base b representations of k and x. Then we define the kth. Walsh 
function by 

walfc(x)=u;r-^^^"^'^"'-'"-'l[o,i)(^), 

where ujb = e^'^^^^. 
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We define translations and dilations of wal^. For j,l G 1: let 

Wj,k,i = b'^^^walk{b'^x - I). 

Note that the support of Wj^kd is given by [bH, Vil + 1)). 

The system {wj^k,i '■ k G No,j, / G Z} is overdetermined in L2(M). Nonetheless, one can 
identify subsets which are complete orthonormal systems in L2(M). 

Let J, /c, / be integers with fc > 0. To each function Wj^k,i there corresponds a tile T given 

by 

T = T,- = [VI, V{1 + 1)) X [b-^k, b-^{k + 1)). 

For the convenience of the reader we prove some results concerning the orthogonality 
and completeness of Walsh functions. See also |15] and the references therein, where these 
results and further information can be found. The following results will be sufficient for our 
purposes here. 

Lemma 1 Let j, j' , k, k' ,1,1' be integers such that k,k' > 0. Then Wj^k,i is orthogonal to 
Wj',k',v if and only ifTj^k,i n = 0- 

Proof. Assume that 

Tj^k,i n Tji^k',v = 0- 

If 

[Vi,V{i + i))n[v'i',v'ii' + 1)) = ^, 

then the functions Wj^k,i and Wji^k',i' have disjoint support and hence are orthonormal. Assume 
now that 

[Vl,h>{l + 1)) n \V'l',V'{l' + 1)) = [max{&'7,6^'r},min{6^(/ + + 1)}) ^ 

and 

[b-^k, b-^{k + 1)) n [b-^'k', b~^\k' + 1)) = 0. (1) 

We consider the case j < j' (the other case can be shown analogously). Then \})~^k, b~^{k+ 
1)) C [b-^'k',b~^'{k' + 1)) and therefore 

max{6^7, V'l'} = Vl 

and 

min{&^(/ + l),V'{l' + 1)} = V{1 + 1). 
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Therefore, using the substitution x = b' {y + I), we obtain 



'Wi^k,i{x)wj'^k',i'{x) dx = fc"-?/^^-? /2 / walk{b~^ X — l)-walk'{b~^' X — I') dx 



-'3,' 

JbH 

= b^j-j'y^ f walfc(?/) W8\k>{yhj-'j' + Ip-^'l - V) dy 
Jo 

= b^^-^'^^^walk'ih>-^'l-l') [ walkiy)walk'iyb^-^')dy 

Jo 

= b^^-^'y'w^h'{h>-^'l-l') [ waU(y)walL,,y./j(i/) dy. (2) 

From (II]) it follows for j < j' we either have b^^'k' > b^^(k + 1) or b^^'{k' + 1) < b^^k. 
In both cases we have k ^ \ k'V^^ J and hence ([2]) yields 0. 

Now assume that Tj^k.i^Ty ,i' 7^ 0- Then the support of Wj^k,i and wy^w^v is not disjoint. 
Using the same arguments as above we arrive at ([2]). Assuming again j < j', we have 
[b-^'k',b-^'{k' + 1)) C [b-^k,b-^{k + 1)) which implies k = [k'b^-^'\. Hence yields 

w,,k,ii^)wf,k',i'{x) dx = b^^-^'^/^wah'ib^-^'l-l') ^ (3) 

and the result follows. □ 
Lemma 2 Let j, k,l & such that k > and b\l. Then 

span {wj^k,h uij^k,i+i, • • • , Wj,k,i+b-i} = span {wj^i^kb,i/b, Wj+i^kb+i,i/b, • • • , Wj+i^kb+b-i,i/b}- 
Proof. Using ^ we obtain for < r < 6 that 

6-1 6-1 

y^(^^j+l,fcfe+r,Z/fe, Wj^k,l+s)L2 Wj,k,l+s{x) = b~^^'^^Wa\r{s/b)Wj^k,l+s{x) = Wj+i^kb+r,l/b{x) 
s=0 s=0 

and 

b-l 6-1 

^{Wj^k,l+r,Wj+i^kb+s,l/b)L2l^j+l,kb+s,l/b{x) = b'^^^^Wals{r/b)Wj+i^kb+s,l/bix) = Wj^k,l+r{x), 
s=0 s=0 

hence the result follows. □ 

The dual of Lemma [2] in terms of the corresponding tiles can be stated in the following 
manner. Let j, k, I be as in Lemma [21 Then: 
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the tiles Tj^k,h Tj^k,i+i, • • • , ^j,fc,«+fe-i cover the same area as 
the tiles Tjj^i^kb,i/b,Tj^i^ki,^i^i/i„ . . . jTj^i^j^h^i^^i^/i,. 

Lemma 3 Let r and t' be two finite sets of tiles such that all pairs of tiles in r are disjoint 
and all pairs of tiles in t' are disjoint. Let W and W be the corresponding sets of Walsh 
functions. Then 

[jT= [j r ^ span = span Vr'. 

Ter T'er' 

Proof. The proof follows by successively using Lemma |2j □ 
Lemma 4 Let J C Z, K C No, and L C Z and set 

r = {T,^ky.jeJ,keKJeL}. 
Assume that the tiles in r are pairwise disjoint. Then the system 

W = {wj^k,i : j e J,k e K,l e L} 
is a complete system in L2(M) if and only if 

y T = M X M([. 

Ter 

Proof. If 

y T = M X M^, 

Ter 

then by applying Lemma [3] one can obtain all tiles of the form TQ^k,i, G No and / G Z as a 
finite linear combination of tiles in r. These tiles also cover M x M.^ , and the corresponding 
set of Walsh functions includes the classical Walsh function system, which is known to be 
complete in L2(M). Since the span stays unchanged when applying Lemma [3] it follows that 
the system 14^ is a complete system in L2(]R). 
Assume that 

y T ^ M X M+. 

Ter 

Then there is a tile T' such that T' C M x M.q \ IJj^g^ T, then the corresponding Walsh 
function is orthogonal to all functions in W and hence W is not complete. □ 
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In the following we consider functions / : — )■ M. In this case we define tensor products 
of the Walsh ftmctions in the following way. Let j.l G and /c G Nq be given by j — 
(ji, . . . ,js), I = ih, Is), and k = {h, ks). Then 

s 

1=1 

The tile corresponding to 'Wj^k,i is given by 

s s 

Tj,k,i = W^jM = n + 1)) X [b-^^k,,b-='{k, + 1))) . 

i=l 1=1 

All the results of this section also hold for the tensor product case, that is, two functions 

'^j,k,i,Wj' i' are orthogonal if and only if the corresponding tiles are disjoint and a system 
of orthogonal functions {wj^k,i '■ G R} is complete in L2(]R*) if and only if the 

corresponding tiles cover (M x M^)*. 

3 Smoothness, convergence behavior, and the decay of 
the Walsh coefficients 

In this section we define classes of integrands of functions / : ^ R. The smoothness of 
the integrands will be controlled by local smoothness parameters and the rate of decay of 
f{x), as the point x tends to infinity (in one or more of its coordinates), is controlled by 
local weight parameters. 

Let s > 1 and let Vr^ be a partition of R* into subintervals J of the form J — [IP I, V -k 
{I + 1)) C where j,l G V. 

We control the rate of decay of the Walsh coefficients of the integrand / using three 
parameters: 

(i) The local smoothness of the integrand /, denoted by for 7^ m C 5* (we assume 
1/2 <au< 1); 

(ii) The rate of decay of the integrand / as |a;|oo 00; for this we use a function 70 : 
V^s R+; 

(iii) The rate of decay of the 'derivative' of the integrand /; for this we use functions 

7„ : Vrs ^ R+ for 7^ m C 5; 

We introduce some necessary restrictions on the functions 7^. 
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Definition 5 Let s > 1 and let V^s be a partition of into suhintervals J of the form 
J = [IP ^l^ip i,{lJ^ 1)) c W, where j, I E . 

Then we call 7 = {'ju)ucs local weight parameters if the weight functions 7^ : Vrs — )■ M+ 
are such that for each u O S we have 

sup 7„(J) < 00. 

We call a = {au)(D^ucs local smoothness parameters if the functions au '■ Vrs — ]R+ are 
such that 

sup au{J) < 1 (ind inf a„(J) > 1/2 

fordl^uCS. 

Note that the assumption infjgpjj^ Q;„(J) > 1/2 is needed for the main results of the paper, 
hence we include it already in Definition O 

We now define a local variation for functions / : R"* — )■ M. Let J = [tP-kl, tP-k{l + l)) C R"* 
for some j,l € Z*. For a subinterval I = Y[i=i[^iyyi) — J with Xi < yi and a function 
/ : R* — 7- R, let the function A{f,I) denote the alternating sum of / at the vertices of / 
where adjacent vertices have opposite signs. (Hence, for instance, for / = ^1^=1 /« have 

Aif,I) = UU{f^{x^)-f^iy^))■) 

Let J = [b^ T*r ^ (Z + 1)) C W for some j,l E V . We define the local generalized 
variation in the sense of Vitali of order l/2<Q;<linJby 



y^Af) = sup 5^ voi(/) 



A(/,/) 



Vol(/)^ 



1/2 



where the supremum is extended over all partitions Vj of J into suhintervals and Vol(/) 
denotes the volume of the subinterval /. (Again, one could include the cases where < a < 
1/2.) 

For a = 1 and if the partial derivatives of / are continuous on J we also have the formula 



d^f 



dxi ■ ■ ■ dxg 




For 7^ u C 5, let V^^j\fu]u) be the local generalized Vitah variation of order 1/2 < 
a < 1 in J„ = -k IP^ -k (/„ + !„)) of the |u|-dimensional function 



fu{,Xu) = / f{x)dxs\u, 
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where Js\u — [b^^^"" * ls\u: b^^""" * {ls\u + l5\u))- For u — $ we define 

\ 1/2 



v;7(/0;0)-(/j/Nrdx) . 



Let Vt^s be a partition of R'* into subintervals of tlie form J = [IP -k l,th> -k [I + 1)). Tlie 
generalized Hardy and Krause variation witli local smoothness a and local weight 7 with 
respect to the partition Vrs is defined by 



,-V2 



A function / for which Va,'y{f) < 00 is said to be of bounded (or finite) variation of order 
a. Further we set 

Hoc,'y = {/ : ^ R : / is continuous and Va^^{f) < 00}. 

Let / : — >■ R be given such that Va-yif) < 00. Let the Walsh coefficient fj k i be given 

by 

fj,k,l = {f,Wj^k,l)L2- 

For j, I e and r e let 

\ 1/2 



2 



'^3,r,l{f) =1 E ■ ■ ■ E l-^.?'*''' 
\fei=[6'-i-ij ks = [b'rs-ii 

where r — (ri, . . . , r^). 

Lemma 6 LetVRs be a partition ofW into subintervals of the form J = [6^ 6^ ★ (j + 1)), 
let ot be local smoothness parameters, and 7 be local weight parameters. Let / : R* — >■ R &e 
given such that V^,-y(/) < 00. 

(i) For any j, I e such that J = [6^' ★ Z, 6^' ★ (/ + 1)) e Vrs and r e \ {0}, where 
Ur — {i & S : ri ^ 0}, we have 

(^3,r,l{f) <{b- l)(""-(-^)-V2)l«r|j,-a..(J)|r|ij,a„,(J)Eie.,Ji^,-Ei^„,jV2^^(J)y^^^ 

(a) For any j,l such that J — [b> -kl^b' -k {I + 1)) e Vr^ and r e Nq, where Ur — {i E 
S : ri ^ 0}, we have 

<Tj,r,i{f) < 6l"-l2l"^l70(J)l-«,^(/). 
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Proof. We prove (i). Let r G \ {0}. We have 
fj,k,i = 



f6^i(«i+l) 



5-(ii+-+i.)/2 



bi' ih+l) 

feJs Is 



/(x)walfc(6 ^ -k X — I) dx 



'[0,1] = 

where -kx = {b~^^xi, . . . , h~^''Xs). We have = {z G 5* : 7^ 0} 7^ and therefore, using 
[2[ Lemma 13.23], we have 

2 



a: 



•,,.,«(/) < (6- i)(2"-i)KI^-i+-+i=5-2"l^li 



< (6- l)('"-^)l"'^lr2°l^l^62aE,e..J-«r^-e..i-^^(J)2y2 (^^^ 



We now show fu) for r = 0. Notice that for fc = we have 



w. 



hence 



<^3,o,i - \fj,o,i\ - 

= 5-(ii+-+i.)/2 

< 5-(ii+-+i.0/2 



/(a;)wj,o,i(a3) da; 



f{x) dx 



[W,ti,fejjr{< + 1)) 



|/(a;)|dx 



1/2 



l/Nrda; 



1/2 



< 5-(ii+-+i.)/2 / Ida; 

\J[bUlfiU(l+l)) J \J[b^*l,bim+l)) 

Now assume r 7^ and let = {a = (oi, . . . , a^) G Nq : < < 6''* for 1 < i < s}. 
For some a G let x G [6^ ^ir Z + V^"^ -k a,lP -k I + IP'"^ i< {a+ 1)) and 

ki=0 ks=0 
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/ f{y)wa\k{b ^ ★ cc - Z)walfe(6-J ★ y - dy 



_ fj-ji js+r i-i hrs^ 



where 



Cr,a^ / f{y)dy. 

J[bJ*l+I>>—^*a,bj-kl+bi—^*{a+l)) 

Let now gf(a;) = for a; ^ [b^ -kl,V -k {I + 1)) and otherwise let 

= I](-l)'"'^i,r-(l.,0^\J,«(a;)- 

Then (? is constant on intervals of the form [IP -kl + -k a^lP -k I + J>>~'^ ★ (a + 1)) and 
therefore 



= / \9{x)\'dx 

J\b3-kl,bi*(l+l)) 



^^y^^ J [b3i,l+b3-'^*a,b3*l+b3-ri,(a+l)) 



— lj-{ji+-+js)+ri+-+rs 



E (-l)'"'cr-(l„,0s\„),(La„/6J,as\„) 



where (Lau/6J, a5\„) is the vector whose ith component is [aj/^J for i e m and Oj otherwise. 
Let 



dr,a = / 1/(2/) I dy. 

'[bi-kl+bi-^*a,bi-kl+bi-^*{a+l)) 
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Then 



a&Ar 



< l)-Ui-^ ^js)+ri+---+rs^\ur.\ 

X V /" Idy 

X / \fiy)\'dy 

J[f)3'^i+6J-'-+{i«r>o)*[a/6J,feJ,t<+bJ'-^+(i"r-.o)*([a/bJ + (l„^,0))) 

< 



5KI4KI / |/(?;)rd?/ 



52|«.l4l«.| /■ |/(l/)rdy 



□ 



We analyze now the behavior of the Walsh coefficients in terms of 01,7. Let J = [b^ 
/,6^'^(Z + 1)). 

(i) Roughly speaking, the parameter controls how fast aj^r,i{,f) decays when the lo- 
cation J (i.e. the support of the Walsh function Wj^k,i) is fixed but the frequency r 
increases. This follows from (i) of Lemma because of the factor note that 
the dependence of ^^^(J) on the location is limited, since au^{J) < 1. 

(ii) The function 70 controls how fast (Tj^r,i{f) decays when the frequency k (or r) is fixed, 
but the location J changes. This is modelled through the behavior of 70,^,; and follows 
from {a) of Lemma |6l 

(iii) Let u C S he fixed and let r = (r„, 0) with r„ G N'"'. Now consider a change of 
location and frequency in the coordinates in u simultaneously. From Lemma [6] (i) it 
follows that aj^r,i{f) decays with 6~""l''li7„( J). Hence if 7n('/) decreases as J moves 
towards infinity, then the diagonal elements, where the frequency and location increase 
simultaneously, decay faster than if just the frequency increases. 

We prove a result concerning the convergence of the Walsh series for functions / with 
Va,'y{f) < 00. The result is analogous to [H Theorem XVI]. 
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Theorem 7 Let cx be local smoothness parameters such that 

inf a„(J) > 1/2 

for all (1} ^ u S and let 7 be local weight parameters. Let f G Ha^-y- Let D C 1?'^ be such 
that {\b' -k I, IP -k {I + 1)) : (j, Z) G D} is a partition ofW and consider the set of tiles 

{Tj,k,i ■■ [b''^~^\ <ki<¥^forl<i<s for some (j, I) E D,r e Wq}, 

which forms a partition of (M x M^)*. Further assume that for all u S we have 

where. J = [b^ l,b^ {I + 1)) . 

Then the Walsh series 

fj,k,iWj,k,i{x) 

converges absolutely and we have 

fix) = Y Yl fj,k,lWj,kAx) 
{j,i)eD ke% 

for all X eW. 

Proof. From Lemma M and the Cauchy-Schwarz inequahty we obtain 

fci = [bn-ij k,=lb^s-i\ 
< (5 _ X)"«r.("')l«r-|^(l/2-a„^(J))|r|i^a„^{J) j,/2^^ {J)Vcy_-y{f). 

By the assumptions of the theorem, the last expression is summable and hence the Walsh 
series is absolutely convergent. 
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Since the Walsh series converges absolutely, its partial sums form a Cauchy sequence. In 
[13] (or see also |21 Appendix A. 3]) it was shown that the sums 

(i,OeDfceni=i{o,-,&'''-i} 

converge to f{x) as ri, . . . , — )■ oo. Hence the convergence of the Walsh series to / follows. 
□ 

We assume throughout the paper that the assumptions of Theorem [7] are satisfied. 



4 Numerical integration 

In this section we study the worst-case error for numerical integration in the unit ball of 



Ha ,y, that is, 



e{Ha,f,Qp,A) = sup 

/eHo,^,yc«,^(/)<i 

Here, the quadrature formula is of the form 



Qp,a(/)- f fix)dx 



n=0 

where Aq, . . . , Aat-i are positive weights and Xq, . . . , x^^i G W are quadrature points. The 
guiding principle for choosing the weights is the idea to have equal weight quadrature rules 
locally on elementary intervals of which yield a small integration error. Let J be an 
elementary interval from a given a partition of M'' into elementary intervals. Let Nj be the 
number of quadrature points in J, then the weight corresponding to those quadrature points 
is given by Vol(J)A^J^. 

We consider tilings of the phase plane which allow us to use Lemma [61 For a location 
fixed by J, Z G Z** we include all frequencies G Ng. Let D C Z^** be such that the intervals 
[V * I, V + 1)) for (j, I) eD form a partition of W. The set 

B = {{j,k,l):{j,l)eD,keni} 

defines a disjoint tiling {Tjk,i ■ (j, k, I) G B}, which covers (R x Mq)'^, that is, 

• Tj^k,i n Tj^^k'x = for all U, k, I), (/, k', I') G B with (j, k, I) ^ {j', k' , I'), and 
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X 



This ensures that the corresponding system 

{w3,k,i ■■ {j,k,l) e B} 



is a complete orthonormal system of L2 
For {j,r,l) e B let 



n=0 -^^^ 



1/2 



Theorem 8 Lei a be local smoothness parameters and 7 be local weight parameters. Let 
Qp^A be a quadrature rule and let B, Sj,r,i, and <Jj,r,i{f ) be defined as above. For r e Nq let 
Ur — {i & S : ri ^ 0} . Then we have 

U,r,l)eB 
{j,l)eD reNg\{0} 



Proof. Let 

(j,k,l)&B 

Consider the quadrature formula 



a; 



Af-l 



<9p,a(/) = ^ A„/(a;„). 



n=0 



We have 



Q 



PAif) - [ fi^) 



~ Yl fj,k,i 
U,k,i)eB 



dx 

'N-l 
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{j,r,l)eB ki = [b'-l-'^i fc^ = [fe'-s-ij 

Using the Cauchy-Schwarz inequality we obtain 



'N-l 

n=0 -^^^ 



f{x) dx 



The second bound follows by using Lemma (i)- 



□ 



In the bound on the integration error in Theorem [8] the only factor which depends on 
the quadrature points P = {xq, . . . , x^-i} and weights A = {Aq, . . . , Aat-i} is Sj^r,i- In the 
following lemma we show that a certain choice of weights will guarantee that many 5-/,o,i are 
zero. 

Lemma 9 Let P = {xq, . . . , a3jv_i} be a set of quadrature points. For (j, I) E D let 

Nj,i = {0 <n < N : x,,e[¥ + 1))}. 

Then we have Uy.oeo = {0, . . . , iV - 1} and Nj^i n N^^i^ = for (j, Z) ^ (j', I'). For 
n G Nj^i let 

A„ = = I I ■ (4) 

Then for all {j, r,l) E B we have: 
(i) If\Nj^i\ >0, t/ien 5,- o,« = 0; 
(ii) If \Nj^\ = 0, then 5j,o,« = 
(zzz) If \Nj^i\ =0 andr eWqX {0}, then 6j^r,i = 0. 

Proof. Since the intervals [IP -kl, IP * (Z + 1)), for (j, I) G -D, form a partition of M'^, the sets 
Nj^i form a partition of {0, . . . , — 1}. 



We have 

/oo 
-oo 

Hence, for r = and \Nj^i\ > we have 



i+-+js)/2 



n=0 



W'^ ifA; = 0, 
otherwise. 



EUh+-+js)/2 



neN. 



\Nj,i 



0. (5) 
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This shows (i). 



If \Nj^i\ = 0, then the sum "^^=0 ^n'Wj,o,i{Xn) = 0. Hence (jS]) imphes 
For r 7^ and \Nj^i\ = we have 



n=0 



E ••■ E 

TV-l 

n,n'=0 fci = [fe'-l-ij fcs = Lfe''''-iJ 

= 0, 

since no point of P lies in the support of Wj^k,!, which imphes [in). 

For the remainder of the paper we assume the weights A are given by (j 
we have for r ^ that 



□ 



In this case 



ir2 



-Jl Js 



A 



= 5.1+-+.. 5^ ■■■ ^ ^ 5^ wMib-'^Xr,-i)eib-Uxn'-i)). (6) 

Using the bound in Theorem [8] and Lemma [9] we obtain the foUowing result. 

Lemma 10 Let a be local smoothness parameters and 7 be local weight parameters. Let 
P = {xq, . . . , xn^i} be a set of quadrature points and A = {Aq, . . . , Aat-i} be given by (jl]). 
Let Sj^r,i be defined as above. Then we have 

e{H^,-y,Qp,A)< Yl b^''+-+^'^/'^,{J) 
U,i)eD,\Nj^i\=o 

The error bound in the lemma above consists of two parts. The first sum arises from the 
fact that the points are only spread over a finite area and hence one obtains a truncation 
error, whereas the second term arises from the approximation error of the integral over the 
region where there are points. 

In the following we define translated and dilated digital nets in arbitrary elementary 
intervals in R^. 
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Definition 11 Let j,l G be given and let P = {xq, . . . ,Xhm_i} be a point set in I = 
[bj i< I, V * (Z + 1)). Let T :[!>' i< I, V {I + 1)) [0, 1]" be given by 



Then we call P a digitally shifted digital {t,m^ s)-net over in I, if the translated and 
dilated point set {T(a;o), . . . ,T{x},m_i)} c [0, 1)'* is a digitally shifted digital net over Z^. 

In the following we analyze the integration error when for each (j, I) E D with Nj^i > 
the points in [b^ icl , IP -k [I + 1)) form a digitally shifted digital net. 

Lemma 12 Let (j, I) E D be given. Let P = {xq, . . . , a;ftm_i} c [b^ I, V {I + 1)) be a 

digitally shifted digital {t, m, s)-net over TL^, in \fp * Z, \P t^t (Z + 1)). Let r G Nq \ {0}. Then 



Proof. Let = T{xn) for < n < 6"*, then Pt = {Uq, . . . ,yf,m_i} is a digitally shifted 
digital {t, m, s)-net over Z^. We have 



T{x) = b-^ irX-l. 



6. 



if < m — t, 

(1 - l/5)KI/25(ii+-+Js)/25(k|i-m+t)/2 ^Jl^l^ >m-t. 




Z) e {b-^ *Xn' - Z)) 



fci = [6n-ij 



fcj, = [6'-l-ij n,n'=0 




Let C 



, Cs G Z^^™ be the generating matrices of Pt and let 

D = {fc G : CiA^i + ■ ■ ■ + Csh = OeZ]^} 



denote the dual net. Then 




Hence it follows that 5j>,« = Oif|r|i<m — t. 

If |r|i > m — t, then, as in the proof of pj, Lemma 7], it follows that 



6l^^i < (1 - 



which implies the result. 



□ 



20 



Lemma 13 Let (j, I) E D be given. Let P = {xq, . . . , x^^^i} C [b^ * I, V -k {I + 1)) he a 
digitally shifted digital {t,m, s)-net over in [b^ -kl^V * (Z + 1)). Let 1/2 < a < 1 and let 
7^ u C 5. Then 

J2 &""'''"'^<5i,K,o),i <{b- i)-W%H/%m--b^n+---+3s)/2^-a{rn-t) -t + \u\\ 
Proof. We have 

J2 ^~"'''"''^i,K,0),« < (l-l/6)l"l/260-i+-+^-=)/26(— +*)/2 6(V2-a)|r„|i 

r„eNl"l ruGNl«l,|r„|i>m-t 

oo 

= (1 - l/6)l^l/25(jl+-+i-)/25(-'"+i)/2 ^ l)il/2-a)l ^ ^ 

;=m-i+l r„eNl"l,|r„|i=« 

< (1 - l/6)l"l/25(j'l+-+J''')/25(-"^+*)/2 ^ 5(l/2-a)/ /^^ + 1^1 ~ 

where the last inequahty follows from [H Lemma 6]. □ 

The following theorem follows from Lemmas [10] and [131 

Theorem 14 Let a 6e /oca/ smoothness parameters and 7 6e /oca/ weight parameters. Let 
a set of quadrature points P = {xq, . . . , x^^i} be given such that for each {j, I) E D with 
\Nj^i\ > the point sets {Xn : n G Nj i} are digitally shifted digital {tj^i, rrij^i, s)-nets over Z^. 
Let A = {Ao, . . . , Aat-i} be given by dl]). Then we have 

ij,l)eD,\Nj^i\=0 

07^«C5(i,OeD,|Af^,i|>O V I I / 

w/icre J = [6^' 7^ /, i< {I + 1)) anc? 

Ca (J) |«| J = (fe - l)(°"(-^)~^/2)|n|^|«|/2^1/2-Q|„|(J)^(a„(J)+l/2)E,g.i>_ 
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The last theorem lends itself to the following strategy. Use randomized digital nets in 
regions J = [IP -kl, V -k{l + 1)) where 7m(J) is 'large'. In regions J where 'juiJ) is 'small' use 
less (or no) quadrature points (or redefine the regions so that the regions themselves become 
larger). Use variance estimators of the variances of the local integration errors and increase 
the number of points where this variance is largest so that the local integration errors in each 
region are of approximately equal size. Hence the values of 'ju{J) do not have to be known, 
instead one adjusts the quadrature points adaptively using local variance estimators. As the 
number of quadrature points increases spread out to cover a larger and larger area. 

A problem that can occur is that, if the number of dimensions is large, there are too many 
subcubes to consider. For example dividing a region into two parts in each coordinate for 
a problem where the dimension s = 100 yields 2^°° ^ 10^^ subcubes. This is infeasible. To 
avoid this problem one can instead only divide the most important coordinates into smaller 
intervals, leaving the majority of the coordinates without any subdivisions. 

Another problem that can occur with this method is that one overlooks an important area 
where no quadrature points are used since the weight functions are not known in practice. 
In this case the number of quadrature points can increase without decreasing the error. 

A further disadvantage of this method is that the one-dimensional projections are not 
optimal, since the points are chosen independently in each subcube. 

If the adaptive subdivision strategy fails, one can use the approach outlined in the fol- 
lowing section (in this case the one-dimensional projections are optimal). 

5 Construction of digital nets 

We now turn to the construction of quadrature points. Let the number of quadrature points 
be b"^, where 6 > 2 is a prime number and m > 1 is an integer. First, one constructs the 
one-dimensional projections for each coordinate and labels the points from to 6*" — 1, as 
illustrated in FigureO Then one uses a digital net in [0, 1)"* with points a^o; • • • ? Xbm_i, where 
Xn = {xn,i, . . . , Xn,s) and maps it to M'^ using the labels for each one-dimensional projection. 
That is, the point Xn,i is replaced by the point on M with label Xn,ib"^ (note that for digital 
nets, Xn,i = kb~^ for some nonnegative integer k). See Figure |5] for an illustration. We 
present the details of this procedure in the following. 

Let 6 > 2 be an integer. For 1 <i < s and 1 < (ij < Aj let Ji^d, = [b^'''''h.,d^, V''"^^ {h,di + 1)) 
for ji^d,, k,d, e Z such that Ji^d, n Ji^dr = for 1 < < < Aj. Let < m^^A, < "^i,A,-i < 
■ ■ ■ < TTii^i < m be integers such that 

Let 2;ft'"i,i_|_...+fc'"i,di-i, . . . , 'i«-i+fe'"i.<^»-i be b'^^-'^i equally spaced points in Jj^^. such 
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a) 



H \ \ h 



c) 



5 2 3 1 4 6 
-> 1 ' t ■ t f- 



d) 



2 ^ a J ¥ 
t ' t ■ t 





/ 






/ 




4 




t ' t ' 1 



Figure 1: An illustration of the procedure in dimension one in base 6 = 2. a) Decide on a 
partitioning o/M; here we consider the intervals [—4, —2), [—2, —1), [—1, 0), [0, 1), [1, 2), [2, 4); 
b) Decide on how many points should be in each interval and then place them equally spaced 
and left-centered; For instance, the interval [—1, 0) contains the points —1 and —1/2; c) Label 
the points, starting with the intervals with the largest number of points and continue with 
intervals with a smaller and smaller number of points; d) Map the points n/2™, < n < 2™ 
to the point with label n. 



that 



2b'"»,i+...+6'"«,di-i+l 



i,di —mi,di 



S,l+...+6"''>di-i+A: 



where for di = 1 we set b"^''^ + ■ ■ ■ + 6™'''**-! = 0. 

Let Ci, . . . ,Cs e Z^^™ be the generating matrices of a digital {t, m, s)-net over Z^. See 
for instance [21 El El [HI [E] for explicit examples. For 1 < i < s and < n < b"^, where 
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n — no + nib + • — h rim-iV^ , let 



/ no \ / ?7„,i,o \ 

\ rim-l J \ r)n,i,m-l ) 



and set 
Define 



rin,i = Vn,i,ob"' ^ + r/n.i,!^™ H ^ r)n,i,m-l- 



2 



and (xji^i, . . . , Xyi^s)- 



Theorem 15 Let P = {xq, . . . , x^n^^i} be constructed as above based on a digital {t, m, s)- 
net over Z^. For d — {di, . . . ,ds) G ni=i{l) ■ ■ ■ ■> ^i} 

s 
i=l 

Then for each Ja — Y[i=i Ji,di> ^ ^ di < Ai for 1 <i < s with 

rrid > t, 

the points P in Ja form a digitally shifted digital {t.rrid, s)-net over Z5 in Jd- In particular, 
for all d e ni=i{l> ■ ■ ■ i^i] ^il^h rud >t it follows that Jd contains at least one point of P. 

Proof. Let 7 = [6^ ★ qr, 6^ ★ (q + 1)) C such that Vol(7) < b>i,d,+-+js,ds+t-ma^ ^^at is, 

(ii,di - Pi) + • • • + {js,ds - Ps) <md-t. 

We need to show that I contains exactly b"^'^~^^^''^^~^'^^ (is.d^-p^) points of P. 

Let £C„ G /, then Xn,i = z^,^^ e [bP'qi,bP^{qi + 1)) for 1 < i < s. Let k = r]n,i - 6""*'^ - 
b"'^-'^i-\ then Zr,^ . = V'Mi^d, + kiV''-'^-'^''''i and 

l)Piq. < h>iMi^di + kib>'^'i-'^''di < F'iqi + 1) 

for 1 < i < s. This implies that 

hPi-iiMi+rr^i.diq. - b'^'Mi^di <ki< bP'-iiA+"'i,di(q. + 1) - 6"**-''* Zi,^. 
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and 



< Vn,i < fcP'-^''.''. {qi + 1) - h"^^Mi^d, + W^'^ + ■■■ + ir^-'^i-^ 
for 1 < i < s. Hence, for 1 < z < s, the digits 

are prescribed and the remaining digits can be chosen freely. Let ai = m — rrii^di + JiA ~ Pii 

0,0 Vn,i,Oi ■ ■ ■ } Ci,ai — 1 Vn,i,ai — 1- 

Let Ci = (q,i, . . . , Ci^rnV for 1 < i < s. 
Then a:;^ G / if and only if 

/ \ / Cl,0 \ 



Cl,ai 
Cs,l 



n 



Cl,ai-1 

Co 



(7) 



The vectors ci,i, . . . , Ci^ai, • • • , Cs,i, . . . , Cg^as are linearly independent since 



ai H h fls = - mi,dj + -pi) =m-md + - Pi) < m - t. 

i=l i=l i=l 

Hence ([7]) has solutions and therefore there are points of P in J^. 

Since 

s s s s 

i=l i=l 



i=l 



i=l 



there are exactly fe'^d-Oi.di-pi) Os.ds-p^) points of P in J^. 

Let Nd = {n : < n < b"^ and a;„ G Jd}- From above we see that Nd is an afiine 
subspace of and hence the set of points in Jd form a digitally shifted digital {t, nid, s)-net 
over Zfc in Jd- Hence the result follows. □ 
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Note that if < t for some d G ni=i{l; • • • , Aj}, it is possible that Jd does not contain 
any point of P. This follows from the fact that ([7]) may not have any solution. 

Further, it is clear from the construction that if the generating matrices Ci, . . . , of the 
digital net are nonsingular, then each one-dimensional projection of the quadrature points 
{cco, . . . , Xbm_i} yields the points zq, ■ ■ ■ , Zb^-i. 

6 Error bound 

We use the construction of Section [5] and Theorem [H] to obtain an upper bound on the 
integration error. 

Theorem 16 Let D = YU=i ^i, where Di C 1? is such that for each I < i < s the collection 
of intervals 

{[bn,V{l + l)):{j,l)eD,} (8) 

forms a partition ofM.. For 1 < i < s choose Aj different intervals from the set ([8]) repre- 
sented by {ji,d,JiA) e Di for 1 < i < s and 1 < di < Ai, where {ji,d,,hA) Ui,drJi,dO 
all di ^ d[. Let Ei = {{ji^^Ji^J I < di < Ai,l < i < s} , E = []•=! Ei, and 

Let a digital {t, m, s)-net over be given. Let t < mi^^i < ■ • • < < m be integers such 
that 

H h b"^''^^ = 6™. 

Let P = { xq, . . . ,Xhm^i} be constructed using the scheme of Section\^ based on the digital 
{t,m,s)-net over Z^. Let F = e E : ma > t}, where j = {ji,d^, ■ . ■ , js,ds) , d = 

{di, . . . ,ds), and 

s 

md = m- ^(m - m^^dj. 

i=l 

Let A = {\q, . . . , \n-i} be given by (jl]). 

Let ct be local smoothness parameters and 7 be local weight parameters. Then we have 

(j,i)eD\F 
<D^uCSij,i)eF 



md-t+\u 
\u\ - 1 
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where J = [IP -k I, V [l + 1)) and 



Cau{J),\ulJ 



(5_ l)("n(J)-l/2)|M|^|«|/2^1/2-Q|„|(J)^K(J)+l/2)E.6„i>_ 



From the proof of Theorem [H] it follows that only subcubes J = [b^ -k l,b^ -k {I + 1)) are 
taken into account where rrid > t. Hence Theorem [16] also applies if one does not use the 
quadrature points which fall into subcubes J where < t. The quality parameter t of 
the original digital (t, m, s)-net used in the construction described in Section [5] is the same 
in each subcube J where > t. Hence, in this sense, if one uses a digital (t, m, s)-net 
and places the first 6'"'* points of this net in each of the relevant subcubes, then the quality 
parameter of the digital nets in each subcube is also t, hence nothing is gained by placing 
the points in each subcube manually instead of using the procedure described in Section [51 



The aim of this section is to show some examples of how Theorem [TBI can be applied in various 
situations. We do not give a comprehensive overview of how digital nets in W should be 
constructed for applications. This depends largely on the applications at hand and is left 
for future work. 



The classical case of integration of functions defined on [0, 1]^ is included in Theorem [16] 
as a special case. It is related to the results on the worst-case setting in [5J. Therein Haar 
functions on [0, 1) were used (which just corresponds to a particular tiling of the Walsh phase 
plane; the decay of the Haar coefficients there was not related to the modulus of continuity). 
Assume we are given a partition Vrs which includes the unit [0, l)**. Let 



for all M C S*. The function space Ha^-y is then related to a Besov space of functions 



Further choose Aj = 1 for 1 < i < s. With these choices one obtains a function space on 
the domain [0, 1]'' and also qMC rules defined on [0, 1]''. The point set P = {xq, . . . , Xbm_i} 
obtained from the construction in Section [S] is the same as the original digital {t,m, s)-net. 



7 Some Examples 



7.1 QMC rules defined on [0, 1]" 



c>^u{J) = ct for all J G V^s and ^ ^ u ^ 



where 1/2 < a < 1 and 




/: [0,1]^^M. 
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In this case we choose D = F = {(0,0)}. Hence D \ F = ^. We obtain the following 
result. 

Corollary 17 Let P be constructed using the scheme of Section\^ based on a digital {t, m, s)- 
net over Z^. With the choice of parameters described in this subsection we have 

where 

C„^|^| = (5_l)("-l/2)H5l/2-a^H/2_ 

7.2 Rational decay of the weight parameters 

Now we consider numerical integration of functions defined on W, where the weights decay 
slowly to as |a:;|oo — )■ oo. 

Let the local smoothness parameters a be constant, that is 

duiJ) = Oi for all J G V^s and for all 7^ m C S', 

where 1/2 < a < 1. 

For J = YlUK h) let 7.(^) = (1 + min(|a,|, |6,|)2-+i/2)-i^ 

s 
1=1 

for 7^ M C S*, and let 

s 

l,{J) = l[{l + mm{\a,\,\b.\r+'/')-\ 

i=l 

Let 6 = 2 and m; = m — 1 — [//2] for 1 < / < 2(m — 2) and m2m-3 = ^2m-2 = 1- Then 

2m-2 

^ 2"^' = 2[1 + 1 + 2 + 2^ + . . . + 2'^-2] = 2™. 
1=1 

For 1 <i < s and di>l let 

f (0,0) forc/i = l, 

. I (0,-1) forrf, = 2, 

KJtA^hA) \ ((^^ _ 3)/2, 1) for di > 3, di odd, 

[ {{di- 4)/2, -2) for c/^ > 4, ti^ even, 
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m — t + 
|m| - 1 



and 

The intervals Jj^i, Jj^2; <^i,3; • • • form a partition of M. Let 

E^ = iUiAAd.) ■■l<di<2m-2,l<i<s} 
and E = Y[i=iEi l^t 

andZ) = ni=iA. 

Let P be the point set constructed according to Section [5] based on a digital (t, m, s)-net 
over Z2. 

Corollary 18 Let P be constructed using the scheme of Section\^ based on a digital {t, m, s)- 
net over Z2 where m > t + 3s and s >2. With the choice of parameters as above we have 

Proof. We have 

/ -r N r 1 for (i,; = 1, 2 



5 "5 

-1 



(l + 2°^»-''0~ for rfi > 3. 
Therefore, for = ^1^=1 "^iA have 

7„(J.)<2-(2"+V2)EU.M, 

and 

70(^d)<2-("+i/2)S-^-''^«. 

Further we have 

We have wij,^^ = m — 1 — [(ij/2] and hence m — rrii^d.. = 1 + \di/2'\ for 1 < di < 2{m — 2). 
Further we have mj,2m-3 = fni,2m-2 = 1- Thus m — mj^^^ = min(l + \di/2],m — 1) for 
1 < di < 2m — 2. From nid > t it follows that 

s s 

t < rud = m — ^^(wi — rrii^di) = m — min(l + \di/2~\ ,m — 1), 

i=l i=l 
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which imphes that 

s s 

-t>^ mill (1 + \di/2'] , m - 1) = ^ min(ji,d^ + 3, m - 1) < m - t, (9) 



m _ 

1=1 1=1 



since m — irii^d, = 1 + \di/2\ = ji^^, + 3. Since m — t < m and 1 + [(ij/2] > 1 it follows that 
ji,di + 3<m— 1 in ([9]). Hence 

F = EE ■.J2J^<'n^-t-3s 

[ i=i 

Since ji > for 1 < i < s, it follows that F is empty if m — t — 3s < 0. 
We have 

^{Haa-: Qp,a) 
For the first sum we have 



oo 



Y 2~''l^'li < 2' 2""' 1 

{i,OeD\F l=m-t-3s+l iGZs,|j|i=« 

+ s - 1 
s - 1 



< 2' Y 2"°' 



i=r?i— t— 3s+l 

^171 — 1 — 23 



where we used [H Lemma 6]. 

Now we consider the second sum. First note that 



s-1 



md = m- + 3) < m - 3s. 



i=l 

Hence we have 
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nczT? V 1^1 / 



We have 



m — t — 2s 
s-l 



s-l \ s 



Thus we obtain 
eiH^-y, Qp,a} 



2 

1/2 



< 2-"("^-*) '2^(3a+2) + 2(1 + 2^/^)^^ 



and hence the result follows. □ 



7.3 Exponential decay of the weight parameters 

Now we consider numerical integration of functions / : — > R where the functions and the 
modulus of continuity decay exponentially fast. 

Let the local smoothness parameters oc be given by 

au{x) = a, 

for all ^ M C 5, where 1/2 < a < 1. 

For J = ni=ik,&i) let 7j(J) = 2-"^i'^(l«'l'l'''l) and 

s 
i=l 

for u C 5'. 

Let b — 2 and mi — m — 1 — [Z/2] for 1 < Z < 2(m — 2) and m2m-3 — "^2m-2 = 1- Then 

2m- 2 

2"^' = 2[1 + 1 + 2 + 2^ + • • • + 2"^-^] = 2"^. 
31 



For 1 < i < s and di G let ji^d, = 

h,di = 

and 



{di — l)/2 for di > l,di odd 
—di/2 if di > 2, even, 



The intervals Jj.i, Jj^25 ^1,3, • • • form a partition of M. Let 

Ei = {UiA^ li,dj :l<di<2m-2,l<t<s} 

and E = n-=i ^i- Let 

Di = {{ji,d,,h,dj ■ di > 1,1 < i < s} 

and D = YIU Di. 

Let P be the point set constructed according to Section [5] based on a digital (t, m, s)-net 
over Z2. 

Corollary 19 Let P be constructed using the scheme of Section\^ based on a digital {t, m, s)- 
net over Z2 where m > t + 2s and s >2. With the choice of parameters as above we have 



e{H„,^, Qpa) < 23-i2-(™-*) (^^_ + 2^(2"+i)2-"(™-*) \ 
Proof. We have 

Therefore, for m C and = H-Li -^i-Ai have 

lu{Jd) = 2-^-i'"^^(l'»'''il'l'^A+i|. 

Further we have 

We have rrii^di = m — 1 — \di/2'] and hence m — rrii^di = + \di/2'] for 1 < di < 2{m — 2) 
and mj^2m-3 = '^i,2m-2 = 1- Thus m — mj^^. = min(l + \di/2'],m — 1) for 1 < di < 2m — 2. 
From rrifi > t it follows that 

s s 

t < = rn — ^^(wi — rrii dj = m — min(l + \di/2~\ , m — 1), 

i=l i=l 
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which imphes that 

s s 

t> J]mm(l+ [rf,/2l,m-l) = J]mm(l + max(|/i,rfJ,|/i,d^ + l|),m-l), (10) 



m _ 

1=1 1=1 



since \di/2'\ = li^^i + 1 ioT di odd and \di/2\ = —h^di for di even. Since m — t < m and 
1 +max(|/j,rfj, \li^d,, + 1|) > 1, it folfows that 1 + max(|/j,dj, {k^d, + 1|) < m - 1 in ([10]). Hence 

e E : ^max(|/i,dj, 1^^,^^ + 1\) < m - t - s 



Therefore we have 



(j,Oer>\F 



For the first sum we have 



{j,l)£D\F l=m-t-s+l ^ ^ ^ ' 



s-\ 



where we used [H Lemma 6]. 

Now we consider the second sum. First note that 



"^d = m — min(l + \di/2~\ , m — 1) < m — 2s. 



i=l 

Hence we have 



{j,i)eF ^ 1^1 - 1 / 



m — t — s 
s-1 , 
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^ 2~a{m-t)2'ias ^ ^ ^ ] |i^| 

as -min(|/i,dj, Ih^a, + 1|) + a max(|Zi,dJ, \li^d, + 1|) < «• 
We have 



i=0 

Thus we obtain 



^{Haai Qp,a) 

^ 23s-l2-{-m-t) I'm — i\ 2s(2a+l)2-a(m-*) — t — S 



s - 1, 

and hence the result follows. □ 



8 Numerical result 

In this section we present some numerical results. We tested the algorithm for the approxi- 
mation of the integral 



/OO POO /»oo / /»oo \ • 

OO J — OO J — OO \ J —OO / 



As underlying digital net we use the first 2™ points of a Sobol sequence. As in the 
examples in Subsection 17.21 and 17. 3[ we set = m — 1 — [//2] for 1 < Z < 2(m — 2) and 
m2m-3 = rn2m-2 = 1- Then 

2m- 2 

^ 2"^' = 2[1 + 1 + 2 + 2^ + . . . + 2"*-2] = 2^. 
1=1 

The partitioning of the real line can be done in different ways. For instance, let 

ai = erfinv(l - 2"') * X for < / < m, 

where erfinv denotes the inverse error function. The number X determines the size of the 
region where the integrand will be estimated. In the numerical experiments below we consider 
two values for X. Then for 1 < / < m — 1 set 

J21-1 = 0,1) 
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and 

J21 = [—o-h —o-i-i)- 

In each interval Ji we define 2™' equally spaced and left-centered points, which we use as the 
one dimensional projections. 

If a point Xn G 11^=1 "^k^ then the corresponding weight A„ is given by 



i=l 



where Ni^^i^^i.^ is the number of points in the interval 11^=1^*5 (Note that if Ni-^,i2,ia < 
2*, then the corresponding weight can be arbitrarily chosen in the interval [0,ni=i('^'i ~ 

A Matlab implementation of this algorithm can be found at 

http: / / quasirandomideas.wordpress.com/20 10/11/11/ qmc-rules-over-rs-matlab-code- 
and- numerical-example / 

We compare the result using the proposed algorithm with using the inverse normal cu- 
mulative distribution function, that is, we write the integral as 

1=1 I I e2v^['=''fi"^(22;-l)+erfinv(2j/-l)+erfinv(22-l)] 

Jo Jo Jo 
We approximate / by 



Q ^ ]_ ^ g20F[erfinv(2a;„-l)+orfinv{2y„-l)+crfinv(22„-l)] 



n=0 



where Xn = {xn, Vn, Zn) G [0, 1), < n < iV are the quadrature points. 

We use the first 2"* points of a Sobol sequence as underlying digital net in both cases. 

The errors and computation times are listed in Table [H There, e(Rs) denotes the error 
using the proposed method, e(invcom) denotes the error using the inverse normal cumulative 
distribution function, t(Rs) denotes the time required for the computation using the pro- 
posed method in seconds, and t(invcom) denotes the time required using the inverse normal 
cumulative distribution function in seconds. 

The speedup in terms of the computation time for the proposed algorithm comes from 
the fact that the computation of the inverse normal cumulative distribution function is 
computationally intensive and the fact that this step is not needed in our method, but is 
needed in the comparison method. In problems where computing the inverse cumulative 
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m 


e(Ks) X = b 


e(Ks) X = 12 


e(mvcom) 


t(Ks) X = b 


t(Ksj A = 12 


t(mvcomj 


13 


U.139UU1 


1.970174 


0.699538 


0.031 


A A OO 

0.032 


A OO A 

0.289 


1 A 


u.zozzyi 


O.OuuiDo 


u.zyuiUD 


U.UoO 


U.UoD 


U.O / 


15 


0.216679 


0.828577 


0.380826 


0.070 


0.071 


1.257 


16 


0.015490 


0.233993 


0.398023 


0.140 


0.140 


2.414 


17 


0.072803 


0.408627 


0.323633 


0.282 


0.283 


4.732 


18 


0.024119 


0.114013 


0.298877 


0.587 


0.587 


9.752 


19 


0.026249 


0.064150 


0.141170 


1.308 


1.315 


19.460 


20 


0.000056 


0.115068 


0.144152 


2.878 


2.946 


37.755 


21 


0.000002 


0.003248 


0.120121 


6.229 


5.970 


75.6635 


22 


0.000199 


0.002157 


0.096750 


12.250 


12.538 


153.900 



distribution function is not a significant factor, tlien one can expect tlie computation times 
to be more similar to eacli otlier. 

The integration error is also significantly lower using the proposed method. There is 
however a dependence on the parameter X. This value determines the size of the sample 
space, which is an important ingredient in the computation. Using a different partitioning 
numbers Gq, .... a„,_i could yield even better results. In high dimensions, say s > 5, it is 
likely to be advisable to use a different partitioning for coordinate indices beyond, say, 5. A 
complete investigation of good choices of partitionings for high dimensional problems is left 
for future work. 
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